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A SOURCE  SOLUTION  FOR  SHORT  CRESTED  WAVES 


by 

R.  C.  MaoCaay 


Berkeley,  California 
April  1954* 


A Qraen's  function  for  the  boundary  value  problem  arising  in  the 
diffraction  of  short-crested  waves  around  o^t-feslas  of  boundad  cross  saotion 
is  prasantad.  The  diffraction  problem  is  formulated  in  a praoisa  way  which 
assures  the  existence  and  uniqueness  of  a solution.  Iha  Graan's  function 
is  so  cons true tad  as  to  make  possible  a representation  of  the  velocity  poten- 
tial at  internal  points  of  the  fluid  in  terms  of  its  values  on  the  obstaole, 
onus  in  general  reducing  the  diffraction  problem  to  th*  solution  of  a 
Friedholm  integral  equation  of  the  second  kind.  IWo  problems  of  interest 
in  the  theory  of  surface  waves,  th«  production  of  wavas  by  a moving  partition, 
and  the  reflection  from  a horizontal  strip  are  studied  by  means  of  the 
Green's  function.  Numerical  results  are  obtained  for  the  first  problem  and 
indications  of  numerical  procedures  given  for  the  seoond.  In  particular,  the 
strip  problem  is  so  fomulated  as  to  make  possible  the  application  of  the 
variational  methods  of  Schwinger. 
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A SOURCE  SOLS  HOW  FOR  SHORT  CEBSUsD  TtkV E3, 

by 

R=  C,  Mao Cssay 


1.  Introduction 

A difficult  but  fundamental  problem  in  the  theory  of  surfaoe  water  wares 
is  the  study  of  diffraction  around  fixed  obstaoles.  Suoh  a situation  arises 
in  the  theory  of  docks,  piles,  breakwaters  and,  as  an  auxiliary  problem,  in 
the  study  of  chip  motlop*  In  three  dimensional  motion  the  problem  has  been  at- 
taoked  by  Fritz  JohnO ) who  suooeeded  in  obtaining  an  integral  equation  formu- 
lation involving  singular  kernels.  Beyond  this  little  has  been  done  for  the 
ease  of  obstaoles  of  finite  extent* 

The  purpose  of  the  present  report  is  to  present  a source  solution  in  the 
oase  of  the  so-oalled  short-orested  wares,  i.e*,  wares  whioh  are  periodic  in 
a direction  perpendicular  to  that  of  propagation.  Suoh  a solution  is  the 
fundamental  tool  in  an  investigation  of  associated  diffraotion  problems  by 
means  of  Frits  John’*  methods,  mysioaiiy,  it  given  the  velvult;  yoUjiti-1  cf  :• 
line  souroe  on  whioh  the  strength  varies  periodically  along  the  line,  ii&th- 
ematioally  it  enables  one  to  express  the  potential  of.  the  interior  of  the  fluid 
in  terms  of  data  on  the  obstacle  surface. 


In  Section  two  we  present  the  function  and  study  its  properties  as  well 
as  indicate  the  general  diffraction  problems  to  which  it  is  applicable*  la  the 
next  two  sections  we  indicate  how  the  solution  may  be  used  in  two  problems,  the 
determination  of  ths  wares  produced  by  a ware-maker,  and  the  refleotion  co- 
efficient for  a rigid  horizontal  barrier  in  the  free  surface.  A few  numerio&l 
results  are  given  for  the  first  problem. 

The  use  of  a souroe  solution,  or  Green's  function,  is  well  known  for 
other  boundary  value  problems.  In  addition  to  yielding  existence  theorems 
whioh  insure  that  the  problem  is  oorreotly  formulated,  they  may  be  used  to  ob- 
tain some  approximate  solutions.  For  example,  it  is  shown  in  Seotion  4 that 
such  a formulation  enables  one  to  invoke  the  variational  procedures  of  Schwinger. 

An  effort  has  been  made  in  this  report  to  present  ideas  whioh  admit  of 
numerioal  computations.  Some  work  of  this  kind  is,  in  faot,  underway  now  and 
will  be  presented  at  a later  time.  , 

2.  The  wean's  Funoti on 


In  order  to  introduou  the  souroe  solution  in  a natural  way  we  will  dis- 
ouss  briefly  the  general  diffraotion  problem  for  short-orested  waves  in  the 
presence  of  oylindrioal  obstaoles. 

Suppose  an  inccmpi-sslbla,  nor.-viseous  fluid  fills  ths  region, -00 < x <+G0? 


* Supersoript  numerals  refer  to  References  at  »rsd  of  report. 
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Z < + C0,  0<y<a  . The  plane  y * 0 is  to  be  a rigid  bottom  and  y = a 

a free  surfaoe.  Ve  assume  the  motion  to  be  generated  initially  by  a short- 
orested  wave  system,  of  small  amplitude,  progressing  in  a positive  x-directiono 
Such  a system  may  be  represented  by  its  velocity  potential. 


<£)l"  = A COSh  y0  y 6 


lit,  x - ft)  COS 

sin 


kZ-R.{*'V'*}  “s  «Z 


(2.1) 


where  A is  a constant  and  / 0 , w0  are  determined  by 

*=  ~1T  = Xo»anh  7o  a,  w02=/0’  -k2  . (2.2) 

It  should  be  noted  that  Equation  (2.2)  puts  an  upper  limit  on  k in  order  that 
the  waves  progress  in  the  x-dircction,  namely  k < K , and  we  shall 

henceforth  assume  mis  to  bo  the 

How  lot  a cylindrical  obstacle  of  bounded  oross  seotior  be  fixed  in  the 
fluid.  %e  resulting  moti""  once  steady  state  is  reached,  will  again  be  time 
periodic  with  frequency  cr  . Moreover,  we  reduce  the  problem  to  one  in 
two- dimens ions  by  assuming  the  Z dependence  to  remain  in  the  form  k Z . 

Under  these  assumptions  we  wriio  for  the  velocity  potential  of  the  ensuing 
motion, 

<£(x,  y,Z,  t)  = Re  (<£(*,  y)<riCTt  ) kZ  (2.3) 

and  there  results,  for  , the  equation 

+ <£yy  * k*  <t>  = 0 in  fluid  (2.4) 


If  we  denote  by  CQ  the  traoe  of  the  obstacle  on  the  x-y  plane,  and  by 
Gf  that  portion  of  y » a exterior  to  C0,  we  find  that  <p  is  also  subject 
to  the  boundary  conditions. 


4>y  - * = o 

on 

CF 

(2.5) 

4>y  =0 

on 

y = 0 

(2.6) 

4>n 

on 

c0 

(2.7) 

I 


s 


where  n denotes  the  normal  to  C0. 

Ic  order  to  fix  the  solution  uniquely  it  is  necessary  to  specify  con- 
ditions at  infinity.  He  would  expect  the  motion  to  oonsist  of  an  incident  plus 
a reflected  and  transmitted  wave,  i.e.,  a oondition  of  the  form 


<f>  - T G ' 


x— +00,  <f>  - <p({)  — as  x^-00.  (2.8) 


It  turns  out  to  he  sufficient  to  demand  only  the  weaker  oondition 

L~cr>{  C {L’V)-'C'/'  (L,y)l2dy+^|^  (~L,y)  - .C*(-L,y)|*  dy}  =0  (2.9) 
for  some  C > 0 , where  i|/  = <p  - <p{,>  represents  the  soattered  wave.  It  can 

be  shown  that  there  is  a unique  (f>  satisfying  (2.4),  (2.5),  (2.6),  (2.7)  and  (2.9) 
and  that  this  solution  doss  not  indeed  satisfy  (2.8). 

Bie  Green’s  function,  G { x,  y,  x',  y'),  whioh  depends  on  a parameter 
point  { *■ 1 i y ' ) is  tc  be  so  constructed  that  it  expresses  values  or  solutions 
at  (x'  , y'  ) in  terms  of  data  on  liie  obstacle  surface  CQ.  It  is  to  satisfy 
(2.5)  on  the  entire  y ■ & plane,  as  well  as  (2.6),  so  as  to  suppress  reference 
to  valuea  of  the  solution  on  y * 0 and  Cy.  in  addition  it  must  have  a singularity 
of  a certain  type  for  (x,  y)  = ( x 1 , y')  in  order  that  the  operation 


/•♦*>  ra 

J J 'C’xx  + Gyy  - k*  G)  $ (x,  y)  dy  dx 

- _ CZJ  'rt 


reproduce  solutions  of  (2.4).  Finally  G must  oarry  the  behavior  of  the  solution 
for  large  I x I and  henoe  must  satisfy  (2.9). 

The  fundamental  singularity  of  Equation  (2.4)  is  found  by  separation  of 
variables  and  is  given  by  Ko  ( k ■>/  x2  + yz  ) where  K0  is  essentially 
the  Hankel  function  of  order  0 and  first  type  with  pure  imaginary  argument.  For 


our  purpose,  the  important  properties  of 


K0  (Z) 


K0  (Z)  = 0 {9  z ) for  large  Z 

Ko  ( Z ) = A (Z)  log  y + B (Z) 
where  A(Z),  B(Z)  are  regular  for  real 


(2.10) 

(2.11) 


and  A ( 0 ) = 


m-t 


:;j  ..  ••  i ' 


It  o&n  be  shown  th&t  tbs  Groan's  funotlcm  as  defined  is  unique  and  we 
prooeed  to  give  suoh  & function.  Consider, 


G(x,y,  x , y ) = 


,i_  r 

w Jo 
c 


= -x 


cosh  y y'  {K  sinh  y(o-y)  — y cosh  y( o-y) 
y {«  cosh  yo  - y sirh  / o} 


I /*•  cosh  yy  { K sinh  y(o-y)-y  cosh  y( o-y') 


{K  cosh  ya  - y sinh  yo} 


cos  w lx  — x;|  dw 


y > y 


cos  w I x - x*  I d w 


y ^ y (2.12) 


2 2 2 

where  Y = w + K and  C is  a eontour  consisting  of  the  positive  x-axis 
except  for  a small  aomi- circle  around  w0  = v/o2  “ k2  (which  is  real  if  k <K). 
This  funotion  was  given  by  Heins ^ in  its  partial  fraction  expansion  form. 

Making  use  of  the  identities. 


• E 

L I 


I ft 

I * 


cosh  y £ {k  sinh  y (a- £ )- y cosh  y jo-  £ -(K-hy)cosh  yj  cosh  y£ 

y {k  cosh  y a —y  sinh  y a}  — y^K  cosh  ya-"y  sinh  ya} 


- v w’ -rk*  b 

Jo  C0S  wxdw=  Kc(k  V **  + b2  ) , 

expression  (2.12)  may  be  transformed  into. 


e-y  iC-t)  ff-yiC+£) 

r jTy  - + , 


(2.13) 


G (x,  y,  x' , y’  ) = 2tt  K0  ( k V (x  - x^T5  (y  + yr)z  )+  j>~:  K0(  k </  U-x')2  (y+y;)s 
/•*  (K  + y)  cosh  vy  cosh  yy'  £ . / . , 

-Jo  — ■ tr 5S>V  - T »r  005  " lx‘x  ldw  • (2_,4) 

le  proceed  to  verify  that.  G has  the  desired  properties.  From  (2,14) 

Gy  =0  on  y = 0 and  from  (2.13)  Gy  ~ KG  = 0 on  y--o. 

Moreover,  G satisfies  (2.4)  sxospt  at  (x'  , y'  ) where  it  na6  a singularity  ox'  the 
prescribed  type.  In  order  to  study  the  behavior  of  large  | x|  it  is  con- 
venient to  make  some  further  transformations  on  (2.14).  We  note  that. 


I 


Res  | !_  (K+  y)  cosh  y y cosh  y y ' e J_  (K ~ y ) e~  >a  cosh  y0  y cosh  y>  y' 

w:»k  l " y ( K cosh  ya-ysinh  yo)  J " w0  sinh  * a ; + 2 £ 

afnh  2 


thus  we  obtain,  by  shifting  the  contour,  C, 


( X,  y,  y')  = 2V  K0  ( k WwP T(y - y')2  ) + ~V  K0(  k /7W7)2  + fy + '/'}*'  ) 


_ ' (K  + y ] e ° cosh  y0  y cosh  y0  y iWnlx-K»| 


w0  sinh  y0  a 


111 


s!  n n 2 r <1 


-Wo 


/■*  CK  4-  y ) cosh  y y COSh  y y • y.,0 1 x - x'l 
o y ( K cosh  y o - y sinh  y a ) 


(2.15) 


where  C'  is  a contour  lying  in  the  first  or  fourth  quadrant  aooording  ns 
x x ^ is  '0  or  <0  . It  follows,  using  (2.10)  that. 


■ (K  + M)g“*u 
w0  sinh  y a 


/i  ^iwslx-x'l  as  ,x|_co  . 


1 + — 

Sinh  2 y0a 


(2.16) 


Now  the  equation  (2.2)  has,  in  addition  to  the  real  roots  ± yQ  , 


an  infinity  of  imaginary  roots  ± i pn  , n = I,  2 j 


Consequently,  although 


we  make  no  use  of  the  faot,  the  above  transformation  could  be  used  to  obtain 


an  expansion  of  the  integral  term  in  a series  of  the  functions. 


In  the  oase  of  infinite  depth,  the  tranef onoation  may  be  carried  somewhat 


further  in  order  to  yield  an  expression  for  G involving  only  real  integrals.  Fbr 


oonvenienoe,  we  ohange  the  coordinate  system  sc  that  y is  measured  vertically 


downward  fremthe  free  surfaos.  Then  (2. 15)beoomes, 


cos  pn 


1 „/OnIX-xM  < , /\  „-Pn,x-x' 

(y  - y ) e , cos  p vy  + y ) s ' n 


G ( x,  y,  x\  y')  = K0  ( k ylx- x')2 +(y-y')2  ) - <?_K(y  *y'>  <?,’*,x-x'1 


’ ® ( y +K ) _ y ( y + v »'  _ i w ( x - x > ) 

3 yfKly)  * 5 


dw  , (2.17 


whore  now  C*  consists  of  the  positive  imaginary  axis  exoopt  for  a small 
semi-oirole  in  the  right  half  plane  about  the  singularity  at  y = 0,  i.e., 
w = i k . Computing  the  real  part  one  is  led  to, 

G (x,  y,  x',  y)  = 2V  K0  ( k y(x-x' )*  + (>•  y' V ) + ~ K0(  k -V/Tx~^+I7+7F ) 

, 2 1 ’<  - k ( y ♦ y /)  i w ( x - x')  • 

Wo  & & 

-UL  /•<Dry»in  y(y  + y')  + K cos  y(y  + y')  - r(x-x')l  x 

7r  jk  [ v (yi  + K 2 ) * JdT  (2*18) 


with  y — -/ r2-k2  and  principal  values  of  the  integrals  are  meant,  and  where 
use  has  been  made  of  (2.13). 


3.  Wave  Generation  by  Moving  Partitions 


As  a first  application  we  consider  the  generation  of  water  wraves  in  a 
basin  by  means  cf  a vertical  partition.  It  is  assumed  that  the  partition  is 
divided  into  seotions  along  its  length,  each  of  which  may  be  moved  with  an  in- 
dependent amplitude  and  phase.  If  the  seotions  are  sufficiently  small,  the 
distribution  of  hcricontal  velocity  of  the  partition  may  be  given  effectively  an 
arbitrary  shape.  We  suppose  the  entire  partition  to  be  sufficiently  long 
that  we  may  negleot  end  effects.  The  problem  is  to  determine  the  motion,  and 
the  thrust  needed  to  obtain  that  motion,  which  should  be  given  the  partition  in 
order  to  produce  ■ prescribed  wave  motion. 

The  solution  is  contained  in  the  following  result  which  is  essentially 
a generalisation  of  the  work:  of  Havelock ^ and  Koanard^. 


(Bf» 


|fcfflfcYiTrv- ' *?' 


*■  . s* 


8. 

Formula  (2.3)  Is  e speoial  case  of  a more  general  theorem  whioh  we 
will  not  develop  hero,  integrals  of  the  above  type  represent  "souroe  layers" 
analogous  to  those  occurring  in  ordinary  potential  theory, and  formula  (3.3)  is 
merely  a statement  of  the  jump  condition  at  suoh  a layer. 

We  are  now  in  a position  to  consider  the  wave  generation  problem. 

To  simplify  the  analysis  we  assume  the  generator  is  of  the  plunger  type  so 


that  f(y)  = B 


a oonstant.  Xe  denote  by 


v l y,  z, 


t ^ 


the  horizontal 


velocity  of  the  partition  and  by 
Noting  that, 

i <3 


V* 


the  surface  elevation  at  infinity. 


V = 


g d 

our  theorem  enables  us  to  consider  several  oases. 
Case  (1)  E real , 


(3.4) 


C 


B cos  cr  t cos  k z 


vs  = q 

V I t_>  ' 


/T  t fin  1/  * 


- ..c  

“ Wo  /|  + ..UL?-\ 

0 V T slnh  2 Yo  0 1 

2 a B/q 


r>  1 s = 


1 1 . 2 71  « 

""0  \ .Inh  9 -v'  n ^ 


co s ( w0  x - cr  t ) cos  kz 
COS  ( w0  s - CT  t ) sin  kz 


slnh  2^0 


Case  (2)  B imaginary. 


v2c  = B sin  cr  t cos  kz  , r)  i,c  — 


:.s  - 
'Iso 


2 cr  B /q 


w0  ( I + 


2 r,° 


slnn  2 pjo 


) 


2 cr  B /q 


w /i  + 2 r«° . \ 

W°'  slnh  2 Tgo  > 


Sin  (w0  x - crt)  cos  kz 


sin  (w0  x - crt ) sin  kz 


Case  (3)  Superposition  of  v,c  and 


v = v,c  + v?s  = B ccs  (kz  - crt ) 


77  = 7)  1 • C 

'GO  'GO 


% 


8,  S 


2 cr  B/q 


w0 


(1  + 


SJfci. ) 

slnh  2 tjo  > 


COS  (w0  x - kz-cr  t) 


(3.5) 

(3.6) 


. 


x * 

^ r 


5-  - y 


1 


(I. 


The.  first  two  cases  represent  situations  in  which  the  sections  of  the 

t 

t 

partition  move  with  the  ^aas  phase  but  varying  amplitudes.  The  resulting 
wave  motion,  sufficiently  far  from  the  partition,  will  form  a snori-oreetwd 
wave  cystera  propage  ting  outward. 

In  the  third  case,  on  the  other  hand,  the  emotions  move  with  the  same 
amplitude  but  differing  phases.  The  resulting  system  is  now  a piano  wave 


-i 


moving  outward  at  an  angle,  tan  — , to  the  z-axis,  with  frequency 


and  wave  lengtu 

^ o 


2 7T 


II 


denotes  t’n6  amplitude  of  the  motion  of  the  partition,  we 
obtain  from  (3.6)  and  (3.6), 


2 K 


Pa  w0 

for  the  amplitude. 


I + 


To  o 


(5.7) 


V CO  * 

• CD 


sinh  2 y0  o 

, of  the  plane  wave  and  in  case  (3).  How  the 


items  we  are  at  liberty  to  specify  sure  cr  and  k , hence  we  take  as  parameters 


2 rr  g 
k 

Wo 


, and  ko  = ror  calculations  the  important  parameters  are  Ka 


Ka  — 

and  = /3  , the  latter  of  which,  as  we  have  seen,  measures  the  direotion 


of  propagation  of  the  plane  wave. 

For  water  of  infinite  depth,  y0  = K , and  thu3 


K o 


2 v 


is  the  ratio 


of  the  depth  to  deep  water  wave  length  of  long-crested  waves  of  frequency  cr 

, (6)  * — * . — * . » Ko 


This  being  the  oase,  we  oan  enter  Wiegel's' 
place  of  d/L0» 


tables  with 


2 v 


taking  the 


A few  results  are  indicated  in  Figure  1.  In  Figure  2 the  corresponding 
oase  of  a flapper  type  generator,  whioh  is  hinged  at  the  bottom,  is  shown.  It  is 
noted  that  the  wave-maker  beoomes  increasingly  less  effioient  as  one  attempts 
to  increase  the  propagation  angle.  In  fact,  as  0 approaches  45° 

the  ratio  in  (3.7)  approaches  0.  If  one  attempts  to  produce  greater  angles  one 
merely  succeeds  In  generating  so-oalled  "edge  waves"  which  vanish  exponentially 
away  from  the  generator.  A similar  situation  involving  a submerged  circular 


waves  from  & strip  of  finite  width  rigidly  fixed  in  the  free  surface.  Math- 
ematically the  problem  is  the  following 

<£xx  + <£yy  - k2  <£  = 0 in  0 < V < a (4.1) 


-e- 

< 

ii 

o 

on 

y - 0 

(4.2) 

6.  - k<£  =0 

• y 

on 

y = a 

1 x 1 > b 

/ » 

V^.O/ 

<Py  = 0 

on 

y = a 

1 x 1 < b 

(4.4) 

<p  - A cosh  y0  y <?lw<,x 

satis 

fies  a 

condition  of  fora  (2-9), 

(4.5) 

He  remark  at  this  stagr  that  the  analysis  of  Frit*  John  may  be  carried 
over  almost  intaot  to  yield  a fairly  general  existence  theorem  for  the  diffrac- 
tion problem  formulated  in  section  two.  This  theorem,  however,  demands  that 
tne  obstacle  surface  bs  perpendicular  to  the  free  surface  at  their  inter- 
section, a condition  which  is  not  satisfied  in  the  present  problem.  Our 
construction  will  in  itself  constitute  such  an  existence  theorem  for  the  present 
case. 

Consider  the  function. 
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9'imi  e.  u 
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= A cosh  "fay1  &1"'  *'  + T,  cosh  -fty  e1"'  *'  &o  X — +00  , (4.7) 


\b  ( y 

r ' • • 


/ / 
w 
i 7 


..  i 

j < 
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Wq  SiTih  J'q  0 
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CC3h#/-  e]w-x'  f*b  <U)e-]"‘Kd> 
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«-  rn  1 
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= A cosh  vy'p1'*0*'  +R  cosh  y0  y'£>'lw°x/  as  x ”00  . (4.8) 
Thu8  i//  ( x y;)  satisfies  (4.5)  with  R representing  a reflected  wave  and 
1*1  oosh  yo  y1  &iw°*  combining  with  the  incident  wavs  to  give  « trans- 

mitted wave.  That  \(/  also  satisfies  (4.1),  (4.2),  and  (4.3)  may  be  verified  by 
direct  differentiation.  In  order  to  mabe  \p  (x;,  y ) a oolution  of 
our  problem  it  remains  to  show  that  f ( x ) may  be  chosen  in  such  a way 

that  (4.4)  is  satisfied. 

Prom  the  defining  equation  (2.12)  for  G we  find. 


o G . / /.  l_  r® sinh  y y' 

d y'  xi  °i  x ) y — 7r  J0  K cosh  y o - y“sTnh  y~o  ' 

On  differentiation  of  the  identity  (2.13)  with  respect  to 


cos-  w ■U^-x4-4w  s (4.9) 
b we  find  that 


(4.9)  m&v  r>a  written  ad. 

\ / v • - » 


d G 


•j4-  (x,  a,  x',  y')  = + -jn-  ||<0  (k/(7-x')2+(a-y')2  )} 


d r 


! /**  ,r  slnh  y y , -y(Q_/)l  , ,, 

ir  Jy  [k  cosh  y o - y sinh  ya  ^ j COS  w|x-x|dw  .(4.10) 


We  now  state  without  proof  the  fbllowing  result,  part  of  which  was 
essentially  proved  in  seotion  three. 

0 ♦ b 

Leama  1 t The  function  <P(x',y')=  l f(  x)  K0 ( k -J  ( x - x')z+(a  - y')2  )dx  , for  f(x) 

b 

continuous  in  -b  < x < +b  , is  a solution  of  <J>X,  x,  + ^yy-k2^  =0  is  y < z , 


! 
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is  continuous  on  y,  = a , - b < x < -r  b when  defined  as  its  limit  value. 


from  y < o and  satisfies, 

"p  d <i>  (x.v)  , , , 

y — o du'1  ' = w ( * 1 ■ or 

y < G 

ft-om  thia  lemma  and  Equation  (4.10)  we  have. 


for  - b < x < b 


/ 1 a <3  y1 


/»+b 

= f (*)  + J f (x)  K ( x,  x')  dx  + iw0  A cosh  y0  c e' 


where 


. . . i,  I /*®  T y slnh  y o , . 1 , / . . 

K(x,x)—  rr  Jn  |_K  cosh  yo  - y sir.hyo  J cos  w * x I d \ 


_ hL.  fa>  cosh  ya  cos  w | x—  x' | dw 
7T  J0  K cosh  yo -y  sinh  yo 


! 
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Setting  . = O gives,  therefore,  an  integral  equation  for  the 

determination  of  f(x) 

Let  us  investigate  the  nature  of  the  Icernel,  K ( x,  x;)  . We  write, 

. /.  I f®  cos  w|»  — x'|dw  I /*®  y cos  w (k-x;)  f tanh  ya-l] 

KU,x  i — - — • — ur*  / , Z ; — " ' dw 

cyo  1 - -K  c:o  (I  - X tanh  ya)  ( !-  £• ) 

where  C*  , C"  consist  of  the  positive  real  axis  exoept  for  eeai-ciroles  in 

the  lower  half  plane  about  the  singularities.  The  first  tern  then  contains 


whatever  difficulties  may  ariso  for  x = x 


Proceeding  as  in  an  earlier 


calculation,  we  find. 


*')  = t r 


' M 1 ~l k " 


t 2 i i 2 

w + k 


d w + 


where  the  dots  indioate  terms  whioh  remain  bounded  as  x — x , Thus, 
aooording  to  Equation  (2.13),  K ( x,  x')  behaves  like  K0  ( k I x-x'l)  near 
x;  = x , and  has  a logarithmic  singularity  there. 
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We  are  thus  confronted  with  the  integral  equation, 

- i w0  A cosh  ■%  a e w°  x = f ( x ' ) +■  / f (x ) K (x , x‘ ) d x , 

b 

in  which  the  kernel  becomes  logarithmically  infinite  as  x -•  x‘  . Ihe 


(4.11) 


Fredholm  theory  may  ba  applied  and  will  yield  the  existence  of  a continuous 
solution  of  (4.11)  provided  there  exist  no  non- trivial  solutions  of  the  as- 
sociated homogeneous  equation.  To  see  that  no  such  solution  exists,  we  suppose 
there  were  one  and  study  its  properties.  Let  f0  (aO  be  the  solution  and 

form,  the  function, 

&(x\y,)=f  f0(  x)  G (x,  a,  x\  y')  d X . 


d G ( x,  a,  x',  a) 


-KGfXiO.x^o)  — 0 in  — b < x < +b  except  at  x = x* 


where  the  difference  becomes  infinite  in  such  a way  that 

f b f(x)  f ' KG(x,  a,  x' . a)l  dx  = f {%') 

J-  k L v y ' J 


Therefore j 


(4.12) 


(It  is  important  to  note  that  this  relation  holds  also  for  the  functions 
\j/  and  f (x‘)  which  are  to  roIit*  our  problem,  a fact  wnioh  makes 

possible  a direct  determination  of  i L{x'.r.)  once  f(x/)  is  known,  without 

recourse  to  another  integration.) 


ft-om  the  expression  (2.14)  for  G,  and  Lemma  (1)  we  conclude  that  ^(x*,  y ' ) 
is  continuous  in  0 < y'  < 0 , ~b  < < b . Moreover  f^x*)  satisfies 

(4.11)  with  the  left  hand  side  set  equal  to  0,  hence, 

iim  0 (*',  y')  _ o, 

>'t°  0 >'  ° ' 
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Finally  ws  make  use  of  another  lemma  which  v*  also  state  without  proof f 
Lemma  Zi  Let  cp  { x,  y)  tea  solution  of  <£xx  + <£  , - k2;£  — 0 in  0<y<a 


satisfying  the  ♦’ollowi.ng  conditions t 

a)  6 y — 0 

b)  <fi  y - K<f>  = O 


on  y = 0 

on  y = a , I x I > b 

on  y = a I x I < b 


c)  y = 0 ony=a  lxl<b 

d)  l — 40  {X  *^x  " 1 C<£{L,y)!2dy  + f I y)- i C </>  (-l-,y)l2  dy  ■=  0 


e)  lim 

« — 0 


/ <£<Pn  dl  = 0 


where 


consists  of  the 


hemispheres,  (x  ±b)2  + y2  = e , y < a 


Then  cfj  ( x , y ) - 0 


in  0 < y < a . 


The  lemma  is  intuitively  obvious,  sines  it  merely  states  that  if  there 
is  no  incident  wave  and  the  strip  is  held  rigid,  no  motion  oan  be  produood  in 
the  fluid.  A proof  may  be  derived  from  an  analogous  Isama  giveB  by  Frits  John. 

The  function  </{,  (x^y*)  enters  the  hypothesis  of  the  lemma,  in  face  will 
be  oontinuou3  at  (±b,a)  j thus  it  is  identically  0 in  0 < y/<0  , and 


by  the  continuity  also  for  y = a.  Sfe  conclude 


f0  (*  > = 0 


in  -b  < x < b 


It  is  interesting  to  note  that  we  may  pass  to  the  limit  as  k — - 0 , 


while  still  retaining  she  same  type  of  behavior  for  large 


I . The  fact 


does  not  seem  obvious  since  the  behavior  of  solutions  of  Laplaoe’s  equation  »t 
large  distances  is  much  different  than  that  of  solutions  of  + (/) yy-k2c£  = 0 

for  »•  > 0 . For  k = 0 we  obtain  for  the  potential  <£(x*,  y')  , 


, , i r*'-  f f®  v.  ' 

*(.',0--*^  m.)> j,  1Kas:;:;i 


where  f ( x ) satisfies. 


— — “T  COS  w ix_x'i  dw  ux  i * iVi 
oh  yd/ 


I 

-iwo  A cosh  wo  ae'*0*'  =f(x')  - «x)[^  K cosh  K?-w  SsTwa  cos  w lx-x'ldw]dx  (4.14)  | 


wo  being  & solution  of 


K cosh  w0  a = w0  sinh  w0  a 
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15. 


If  the  depth  is  infinite  Cha  kernel  anj  bo  evaluated  explicitly  and 


yields. 


P 


-i  wn  be'*0*  = f ( x;)  - -^r  J ( (x)  ^ 7T i e 'IS'*"X ' + cos  K I x-  x I Ci  V K ! x -x;  i)  | 

+ sinKlx-Vl  ( Si  ( K I x -x; !}-•?-)  j dx  (-=--/  j 

vhsrs  Ci,  Si  are  the  integral-oosine  and  :.ntegral-sine  functions  respectively.  j 

. I 

The  appearance  of  the  term  Ci  ( K !x~x  I)  confirms  our  earlier  remark  that  the 
kernel  baoomoe  logarithmically  infinite  as  x —■  x j 

X 

Ihe  existence  of  a solution  having  been  established,  one  may  turn  to 
the  question  of  obtaining  numerical  results.  An  equation  of  the  type  (4.11) 
can  be  solved  in  a straightforward  way  by  replacing  f (x)  by  a polynomial 
approximation  and  thus  obtaining. a sys  m of  n linear  equations  for  the 

determination  of  the  value  of  f ( x ) at  n points.  The  numerical  work  ia= 
volves  numerical  computation  of  integrals  involving  K(x,  x)  and  powers  of  x . 

Ihe  singularity  gives  no  difficulty  as  it  may  be  subtracted  off  and  the  in- 
tegrals involving  it  carried  out  explicitly.  Some  work  of  this  type  is  under- 
way and  it  is  hoped  may  be  presented  in  the  future. 

A more  elegant,  although  possibly  less  useful  method  is  the  application 
of  the  variational  methods  of  Sohwinger^7) , .which  wer®  developed  to  treat 
aperture  diffraction  problems.  If  Equation  (4.11)  is  multiplied  by  Hx) 
and  integrated  over  the  strip  and  use  is  made  of  (4.8)  we  find. 


J Jgal-  /v  / , 

r |K-f  y0)£~  r«u  V ' 
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5l»lll  C.  /Q 


A 

"ft 


f**  f(x')z  dx'  4-  f**  /_;  f(x)  f(x')  K(x,x;)dx  dx' 

n?fu'>  dx'  f 


(4.16) 


where  r is  essentially  the  reflection  ooeffioient.  Now  the  quantity  on  the 
rif£i+  hand  side  of  (4.16)  is  stationary  with  raspeot  to  first  order  variations 
of  fix)  about  its  correct  value  ao  determined  from  the  integral, equation  (4.11). 
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It  might  be  hoped,  therefore,  that  the  reflection  coefficient  would  be  relatively 
Insensitive  to  errors  in  f ( x ) . 

More  explicit  use  of  the  etatl o««ry  eha?"''+er  of  -jr - may  be  made  as 
follows.  Suppose  wo  expand  * (x)  in  (~b,  b)  in  a Courier  series*, 

f (x)  = Y.  <•*•'  6 x 

-CO 

If  we  substitute  this  expression  in  (4.11)  we  obtain. 


(2b  Y.  a.,,  +XI  an  Cmr,)  r ^ ( £ an  Bn  )J 


i t 


where 


♦ D 

1:1  = JJ  6 Mx  6 ""x  K (x,x  ) dx  dx  , Bn  - J 


- in  x , ! w,  x . 

6 S’  dx  . 


Differentiating  with  respeot  to 


and  noting  tiiai.  X is  stationary. 


( 2 a-m"t~Xan^mn)r  — 2 Bm  ( 2^  ^n) 

-CO  1 ' -w 

Defining  constants,  ' Dp  , by 

Qn=-r  Dr>(  T an  Bn)  = ~Q  Dp 

r 1 *"  ' TR  > 1 


m =0,  ± I , ± 2, . . . . 


we  obtain  finally 


P — C,  ±1,  ±2,  . . . , 


(4.17) 


r — X Dp  Bp 

-CD 


(4.18) 


n + > C r>  = B 

wrrr,  w n Wn 

-r 


m = 0,  ±1,  ±2,  . . - 
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We  have  thus  transformed  the  problem  into  solving  an  infinite  system 
of  linear  equations,  (4.19).  Onoe  the  Dn  ’s  are  found  ,the  reflection  co- 

efficient oan  be  computed  from  (4.18).  In  theory  one  oould  also  obtain  f (x)  . 

^ r -»  ' 

The  choice  of  functions  nv  | for  the  expansion  of  f (x)  is  not  es< 

sential  to  the  method.  X h suggests  itself,  however,  since  for  oertain  values 
of  w0  , the  set  i Bn  \ reduces  to  & single  term. 
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Substituting  the  series  for  fix)  with  up  replaced  hy  •~iV  Dp  in 
the  integral  equation  fixes  the  constant  G0  and  thus  determines  On 
For  the  i,ass  of  diffraction  of  sound  waves  through  a circular  aperture,  the 
scheme  has  proved  quite  suooessful,  and  it  was  found  that  only  a few  terms  of 
the  infinite  series  suffice  to  give  quite  accurate  approximations, 
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